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Abstract 

We develop a theoretical model for transverse dynamics of a single electron spin interacting with 
a nuclear spin bath. The approach allows a simple diagrammatic representation and analytical 
expressions of different nuclear spin excitation processes contributing to electron spin decoherence 
and dynamical phase fluctuations. It accounts for nuclear spin dynamics beyond conventional pair 
correlation models. As an illustration of the theory, we evaluated the coherence dynamics of a P 
donor electron spin in a Si crystal. 
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I. INTRODUCTION 



An electron-nuclear spin coupling substantially affects electron spin dynamics in solids. 
This phenomenon is broadly utilized in EPR probing of material structures.— However, for 
novel technological applications of electron spins^i^ it produces a major obstacle even at low 
temperatures when effects of phonons are suppressed.— i^*^^ In particular, it is crucial for 
a quantum coherence that is a key issue of a quantum computation.^ Entanglement of an 
electron spin with a nuclear spin bath results in an irreversible loss of coherence. Unlike spin 
relaxation,— decoherence process is not suppressed even in strong magnetic fields. There are 
several methods to reduce effects of spin bath, such as isotope purification,^ dynamical po- 
larization of nuclear spins or dynamical decoupling of an electron spin evolution. However, 
it is not clear if experimentally achievable values of isotope purification, nuclear polarization 
or precision of electron spin control are sufficient to suppress effects of spin bath beyond 
the required threshold limit. ^■^ Moreover, for some techno logically- important materials these 
methods may be inapplicable. For example, isotope purification cannot be used in GaAs 
nanostructures, because all stable isotopes of it have non-zero nuclear spins. In this connex- 
ion theoretical models of electron-nuclear spin dynamics can provide better understanding of 
electron spin decoherence processes and also help in estimating the effectiveness of coherence 
control schemes. 

In this work we investigate dynamics of a localized electron spin interacting with a nu- 
clear spin bath at a low temperature regime. We chose a system with a long spin relax- 
ation time that is of interest to quantum computing. Until recently, this problem has been 
studied using stochastic models of spectral diffusion .i^^J^iii Their results were verified by 
numerous experiments carried out on macroscopic samples. Recently, emphasis of exper- 
imental and theoretical studies has been shifted to dynamics of single quantum systems, 
where stochastic models are inappropriate. Several analytical and numerical approaches 
based on quantum dynamics have been developed and used to investigate different aspects 
of the problem . Among the issues addressed in these studies 
are spin relaxation at low external fields,— effects of nuclear spin polarization on electron 
spin dynamics,— dynamical control for spin decoherence,— contributions of high order nu- 
clear spin correlations into an electron spin echo,— etc. However, many questions are still 
open. How do stochastic and dynamical models relate to each other? How does one char- 
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acterize short time qubit evolutionY-Stt What are the reversible and irreversible parts of spin 
dynamics?^^ What are the relative contributions of different correlated nuclear spin clusters 
in electron spin dynamics? How do nuclear spin correlations grow in time?— 

Here, we demonstrate that diagram techniques developed previously in studies of Heisen- 
berg ferromagnet o^^i^^i^^i^^i^'''i^^i^^i^°i'^'^ can be applied to evaluate effects of a nuclear spin 
bath on a single electron spin in a high field regime. Our theoretical approach provides 
a transparent representation of different nuclear spin dynamical processes contributing to 
an electron spin evolution. It naturally accounts for nuclear spin excitations beyond pair 
correlation models. We show that transverse evolution of an electron spin can be factorized 
to a precession in a nuclear Overhauser field and more complex dynamics due to electron- 
nuclear spin entanglement.— A conventional Hahn echo experiment cancels the phase due 
to precession in the nuclear field and also suppresses entanglement with the nuclear bath. 
As an illustration, we consider dynamics of an electron spin localized at phosphorous donor 
impurity in a Si crystal. We estimate contributions of 2, 3 and 4 nuclear spin excitations to 
electron spin decoherence and discuss effects of Hahn echo on spin decoherence. 

The structure of the paper is as follows. In the next section we describe a Hamiltonian 
and discuss assumptions used. Section ITTTl is devoted to diagrammatic representation of the 
decoherence process. In Section |TV| we discuss the approach and consider an example of 
P donor electron spin in a Si crystal. Section |3 gives the conclusion. In Appendix El we 
provide the spin diagrammatic rules, in Appendix^ we discuss some specific properties of 
the linked cluster expansion for spin systems and in Appendix O we give explicit analytical 
expressions of some high order nuclear spin contributions to electron spin dynamics. 

II. MODEL 

We consider the spin of a single electron localized in a quantum dot or bounded by a 
donor impurity. We assume that only one type of nuclear spins with / = 1/2 is presented, 
though this assumption can be relaxed within the approach used. In a strong external 
magnetic field the Hamiltonian for a single electron spin coupled by the contact hyperfine 
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interaction-^ with a system of nuclei can be written as 

H = co^S^ - ^/ E + E + E B^^^^i + + Btfl^Ii}. (1) 

A similar Hamiltonian has been used in previous studies of the spectral diffusion problem . 
Here, we briefly describe the notations and assumptions. The first two terms in Eq. ^ 
account for electron and nuclear Zeeman energy level splittings in an external field, H, 
with the Larmor frequencies uje = (7*/3H/h and uj = 7H respectively. The 2;-axis is chosen 
along the magnetic field. The third and fourth terms originate from the contact hyperfine 
interaction. In a strong magnetic field direct electron-nuclear flip-flop transitions, {S^I^^ + 
S^I^), are forbidden by the energy conservation law. Therefore, beside a small visibility 
loss,— this part of the contact interaction contributes to the effective coupling between 
nuclear spins only.— The coupling coefficients are A^^ = (8/3)7r(7e/37|\l/(-Rj)p and BfJ = 
A^Ay /2uje, where Qe = 2 and g* are free and effective electron g-factors, (3 is a. Bohr 
magneton, 7 is a nuclear 7-factor and "^{Ri) is an electron wave function at a position of 
i-th nuclear spin. The last term in Eq. represents the secular part of the nuclear spin 
dipole-dipole interaction. For the electron nuclear hyperfine interaction, we consider here 
only the isotropic contact part. Effects of the dipolar e-n hyperfine interaction, resulting in 
spin echo envelope modulation ,^^i^^i^^i^^ has been discussed elsewhere.— The Hamiltonian (^J 
is diagonal in the electron spin. Therefore, the nuclear spin bath affects transverse electron 
spin dynamics only. 

The initial state of the electron spin plus the system of nuclear spins is described by 
the density matrix p(0) at the time moment, t = 0, when the electron spin state has been 
prepared. Two assumptions are applied to p(0). First, we use the standard approximation 
of a factorized system and bath^^^ p(0) = (8> p°. The electron spin is initially prepared 
in the pure state, {1/ \/2){\+) + |— )), e.g., by a n/2 pulse. The second assumption is 
that the nuclear spin system is in a pure state that is an eigenstate of If operator. 
With the latter statement we neglect by nuclear spin-spin correlations at t < 0. Influence 
of different initial states of nuclear bath on electron spin evolution has been discussed in 
Refs. y and [2^ though the lack of nuclear spin-spin interaction in these papers may affect 
their conclusions. It has been argued that the pure spin state utilized here can be useful 
for quantum computation purposes because it does not destroy electron spin coherence at 
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short timescales. After statistical averaging over possible initial configurations the nuclear 
spin density matrix is = 'Yl,Pn\n){n\, where p„ is a statistical weight of a given nuclear 
configuration \n) = \ ttitii •••)• 

The evolution of the up-down component (H — ) of an electron spin density matrix can be 
written as^^ 

p^_(t) = p° _e-^"'=*Tr„{e-'(^o+^+)Vne'^"^"+^'^*}, (2) 

where 

^o = (l/2)E.^f/f, 

y± = Kid ± Hf , 

Here, we used the fact that the Hamiltonian ^ commutes with 5*^ and projected it to the 
electron spin-up (+) and spin-down (— ) subspaces. The projected operators are written as 
two terms, Hq, which is a sum of single spin operators, and V±, which describes spin-spin 
interactions. A contribution of the nuclear Zeeman splitting, the second term in Eq. (Q), 
is cancelled, because it commutes with the rest of the Hamiltonian. If the nuclear Larmor 
frequency, cu/, varies from a site to site due to inhomogeneity of the external field or other 
factors then its fluctuating part should be included into Hq. Below, we will keep the su- 
perscript indexes, (hf) and (dd), in the coupling constants only if it is not clear from the 
context what type of interaction is used. 
Using the relations 

^-i(Ho+V+)t ^ V+{t'-t)dt'jg-iHot^ 

where T{...} is a time ordering operator, and our assumption on the initial nuclear spin 
density matrix, Eq. is transformed to 
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where Un = {n\2HQ\n) is a contribution of a nuclear Overhauser field to the electron spin 
precession frequency and V±{t) is V± in an interaction picture defined by Hq. For a single 
shot measurement of a single electron spin the nuclear bath contributes to the shift of the 
electron precession frequency, cj^, determined by the initial configuration, and complicated 
dynamics due to coupling between nuclear spins that is described by the bracket {n\...\n). 
The weight factor, p„, corresponds to a statistical averaging over an ensemble of electron 
spins or repeated measurements. For an ensemble measurement, Eq. (jSj) describes a free 
induction decay,— where the transverse magnetic moment decays due to an inhomogeneous 
distribution of spin precession frequencies and also due to spectral diffusion in the presence 
of nuclear spin environment. Eq. ^ can be viewed as an exact formal solution for the 
electron spin dynamics. In following sections we will evaluate the term in it between the 
angular brackets using the Linked Cluster Expansion (LCE) formalism.— 

In Eq. (0) the product of two exponential operators can be transformed to a single expo- 
nential form, T{e~^^-^ ^'■^ }, by shifting the time variable t' —t — > t' in the second exponent. 
One can see that in this case the potential, V{t'), is continuous for the hyperfine- mediated 
interaction given by Eq. 0, while for the dipole-dipole nuclear spin-spin interactions it 
changes sign at t' = 0. To avoid operations with discontinuous potentials we will use a 
two-exponential form of a bracket in Eq. (0) keeping in mind that it can be transformed to 
a single exponent. 

Schematically, the bracket in Eq. (jSJ can be shown as a two branch propagation. Fig. 
The system propagates on a first branch with the Hamiltonian — t) and after that 

on a second branch with —V-{t'). This is similar to the non-equilibrim Green's function 
approach . ^^1^^ 

The factorization of nuclear spin induced dynamics to the phase factor acquired in the 
nuclear Overhauser field and complex dynamics due to entanglement between bath modes 
is consistent with EPR experiments. For example, see Ref. |53|, where an electron spin 
resonance line broadening can be resolved to inhomogeneous and homogeneous parts. 

The above procedure also can be applied to evaluate electron spin dynamics in experi- 
ments where the electron spin is flipped by magnetic pulses. For example, for a Hahn spin 
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echo ( tt/2 — t ~ 7T — t— echo) the evolution equation can be written as 
p+42t) = -p° + ^p„(n|T{e^/o v'+(*-Orf*'}r{e^/o ^-(0^*'}T{e-^/o ^+(*'-*)'^*'}T{e-^ 

n 

(6) 

We emphasize that the phase factor due to precession in the external field plus the nu- 
clear field that appears in Eq. (0) is cancelled for the echo. Moreover, the bracket {n\...\n) 
describing electron-nuclear spin entanglement is different from the one in Eq. These fea- 
tures correspond respectively to elimination of inhomogeneous broadening and suppression 
of spectral diffusion in ensemble measurements.— 

III. LCE AND DECOHERENCE 

Using LCE, we can write the expectation value of a time ordered exponent aa^^ 

(n|T{e/o ^(*')'^*'}|n) = e<^^>+<^^>+<^^>+-, (7) 

where V{t') is an interaction and (Vk) is a contribution of linked diagrams^ only to the 
integral 

[ dt, [ dt2... [ dh{n\T{V{h)V{t2)...V{tk)}\n). (8) 
Jo Jo Jo 

The coefficient, 1/k, that appears in LCE, is included in the (Vk) terms. This expansion 
provides a convenient exponential form to describe the dynamical processes. Moreover, 
each perturbation term (V^) in it corresponds to an infinite sum of terms in a conventional 
perturbation theory. 

The bracket in the expression (jH)) can be evaluated by a diagrammatic technique. The 
term (V^) in Eq. ((7j) is of the A;-th order in the interaction. It describes the collective dynamics 
of a cluster containing up to k spins. Because the proof of LCE can be given based on 
combinatorics, the expansion procedure should be applicable to potentials, discontinuous 
in time, or to products of several evolution operators given in Eqs. (jS}-®- 

Diagrammatic rules for spins are not so transparent as for fermions or bosons, because 
commutation brackets of spin operators do not yield c-numbers. Many papers addressed this 
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issue . In our derivations we use the technique described in Refs. l37ll38lJ41l 
with modifications accounting for specifics of the problem. A brief summary of this technique 
and the used diagrammatic representations is given in Appendix [X] 

The LCE expansion of the scattering matrix ((7j) can be described by the same series of 
diagrams as a free energy in Matsubara formahsm given in Ref. |38l The two-exponential 
representation of the scattering matrix in Eq. (jSJ does not change the structure of the 
diagram series, but affects the spin propagators only. In Figs. [21E] we show the sets of 
diagrams corresponding to the bracket in Eqs. ()5|6p up to the fourth order in the nuclear 
spin-spin interaction. 

We first discuss the scenario where the nuclear dynamics starts from a pure initial states. 
Ensemble results are obtained by taking statistical average of all possible initial configura- 
tions. Dynamics of each distinct configuration of spins in a cluster, in general, should be 
depicted by a different diagram specifying an initial spin configuration and the time arrow as 
it is shown in Fig. |2l However, in many cases analytical expressions for these configurations 
can be transformed to each other, as is discussed in Appendix ^ To reduce the number of 
diagrams, in the figures we omit the configuration dependence, assuming that one diagram 
represents all possible configurations. In analytical and numerical evaluations we calculate 
contributions of distinct spin clusters. We also drop the time arrow for simplicity. 

To evaluate diagrams given in Figs. EHUwe use a 2 x 2 matrix (matrix elements indexing 
branches. Fig. ^ Green's function at an i-th site 

where LOi = /2 and r = ti — ^2- The total evolution time, t, appears in the off-diagonal 
elements of K'(r). Matrix elements of the Green's function (jU)) possess a simple physical 
meaning. The propagator starts at time ti on a branch denoted by a row index and ends at 
time t2 on a branch denoted by a column index. To account for the two-branch propagation 
the spin-spin coupling coefficients, Afj^, Bf^ should be multiplied by a cr^ Pauli matrix (Afj^ = 
Af/a^, Bf/ = Bf^a,) and B}^ should be muhiplied by a 2 x 2 unity matrix (Bg = B}^1). 
In the analytical expressions for the diagrams one has to sum over repeating matrix indexes 
in addition to integration over the time variables. 

The first order linked cluster diagram corresponds to (n| Afj^Jf J||n) (not shown in Figs. El- 
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EJ. However, its contribution vanishes because the sum over the branches is equal to 
Tr{Afj^} = 0. The first non-zero contribution to the decoherence process is due to the 
second order diagrams which represents a nuclear pair spin flips, Fig. |21 The corresponding 
analytical expressions are 

{v^') = - E + ^^}' (10) 

i=Tj=i 

for the hyperfine-mediated interaction only, and 

(vi^) = - E + — ^-^}, (11) 

for the dipole-dipole interaction only, where we define This result is consistent 



with Ref. |25l The real parts of Eqs. pOlll|) contribute to electron spin decoherence while 
the imaginary parts renormalize electron spin precession frequency and, hence, produce 
phase fluctuations. The (V2) diagram with a cross-term, B^^B'^'^, contribution of the dipole- 
dipole and hyperfine-mediated interactions is zero because the spin propagators over different 
branches cancel each other. Therefore, within the second order the contributions of these 
two mechanisms to electron spin dynamics are completely separable. 

We emphasize that in comparison with a conventional perturbation expansion LCE con- 
verges faster. Each order in LCE corresponds to an infinite sum of terms. For example, the 
second order correction in LCE is a partial sum shown in Fig. |S| It includes a series of even 
orders of a perturbation expansion. 

The third order diagrams in Fig. 01 can be divided into two groups. First group includes 
diagrams Efa), Efb) and corresponds to (V2) diagram with attached JfJ| interaction lines. 
Diagrams of this type can be accounted for by renormalization of Green's functions with the 
standard equation. Fig. El where renormalized Green's function (bold line) is 



K^(ti,t2 



e^^"^^(5.i^(-r) - 5,^9{t)) 



[12) 



T = ti—t2, 9 = ti + t2, uj[ = uJi + Aui and Auji = J2j ^jii^j)- This modifies the dipole-dipole 
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pair flip-flop term, Eq. (fTTj). as 



=Tj=i 



(1 -e*^"*)( 




}, (13) 




+ 



1^111-^22 



where uu = ujij + AtUj — A^jj and ^722 = ^ij — AcUj + AtUj. If we assume that Acj, is inde- 
pendent on the site, then Eq. p3|) transforms back to Eq. Therefore, the contribution 
of these renormahzation terms is reduced if the initial polarization of the nuclear spin bath 
is homogeneous. A similar modification of the pair flip-flop term can be due to an inhomo- 
geneous distribution of nuclear Larmor frequencies, uJi in Eq. (Q). The diagram Fig. Efa) 
contributes to renormahzation of a pair dynamics by interaction. By direct evaluation 
one can show that for the dipole-dipole interaction it is also cancelled . 

The second group of diagrams. Fig. Ofc), is due to three-spin flip-flop processes. It 
corresponds to a ring propagation of a spin excitation. For the dipole-dipole interaction 
only, due to symmetry of interaction terms, the clockwise propagating excitation cancels the 
counterclockwise excitation. It is an analog of the Furry's theorem.-^"^ An analytical form of 
the hyperfine-mediated contribution can be written as 



for clusters {i =t, j =i, k =i}, where permutation of j and k spins is included already. For 
=i) J =T) k =]} clusters, one should change signs at frequencies cjjj fc. 
In the third order linked diagrams, an effect of cross-terms including both hyperfine- 
mediated and dipole-dipole interactions appears. For example, a ring diagram. Fig. ^c) and 
diagram Fig. Efa) with two dipole-dipole interaction lines and one hyperfine-mediated line 
have non-zero contributions. Because most of the third order diagrams give zero contribution 
to spin decoherence, the fourth order corrections should be evaluated. 

The set of the fourth order diagrams is given in Fig. |3] The only restriction on vertices 
in this diagrams is that two spins i and j coupled by an interaction line should be different 
{i 7^ j). Therefore, the same diagram in Fig. |3] can correspond to different number of spins 
in a cluster. For example, the ring diagram. Fig. ^f), can describe excitations of two, three 
and four spins. We call them two, three and four-spin ring diagrams respectively. 




(14) 
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In the fourth order in addition to different types of renormahzations of lower order dia- 
grams, Fig.|31^a-e), and a ring diagram, Fig. |3ff) , we have contribution of locked diagrams,— 
Fig. |31^g,h). By locked diagrams we mean diagrams containing vertices with two incoming 
propagators. This group compensates overlapping of spin pair excitations and also restricts 
excitation within a spin space P = ±1/2. The three-spin diagram. Fig. IH^g) modifies a 
double excitation with one common spin, that appear in expansion of the exponent of V2 
term, and also correct the fourth order ring diagram with repeating indexes, see Fig. [7| The 
two-spin diagrams Fig. EJ^g) and Fig. Efh) plays a similar role for pairs with two common 
spins. Fig. |H1 In more details we discuss this in Appendix |Bl 

For spin 1 = 1/2 two-spin diagram is compensated completely by the locked diagrams, 
see Fig. IHb). The tree spin rings describe dynamics ofz=tj=J.fc=|orz=J,j=|fc=| 
clusters. After correction by the locked diagram. Fig. El^b), it corresponds to propagation 
of a spin excitation as i — j — > /c j i. An analytical form of this term together with 
locked diagrams. Fig. Ill^g,h) are given in Appendix O A four spin ring diagram correspond 
to dynamics of three distinct spin clusters tiii; Titi and ttii- 

As a result, up to the fourth order in nuclear spin-spin interactions we write equation for 
the electron spin coherence as 

where the index n denote an initial configuration of the nuclear bath, the (V2)„ term (diagram 
on Fig. 121) is given by Eqs. pOlllj) . the (Vs)^ term (diagram on Fig. IHfc)) with expression 
given by Eq. (fT^ and the fourth order contribution of nuclear spin dynamics schematically 
shown in Fig. 121 Analytical expressions for some terms of (V4)„ are in Appendix O 

The diagram series can be extended to higher orders. At each order there should be a 
group of diagrams renormalizing lower orders with PP terms, a group of locked diagrams 
that compensate overlapping of spin excitations in lower order clusters, and also a group of 
ring diagrams. 

For systems with a high concentration of nuclear spins, contributions of different dia- 
grams to a (Vfc) term can be estimated based on 1/Z expansion^, where Z is an effective 
number of interacting spins. For example, if the last term in Fig. 121 is 0(1), then the first 
term is 0(1/ Z), second term is 0{1/Z) or 0(1/ Z"^) depending on whether the diagram cor- 
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responds to a three or two-spin cluster, and the third one is 0(1/Z^). This follows directly 
from counting of a number of summands (different spin configurations) in the analytical 
expressions, Appendix O If the effective number of spins interacting with a given one is 
large then (V4) term can be approximated by the last diagram in Fig. IHlonly. 

IV. DISCUSSION AND EXAMPLE 

The equation for a free evolution of a single electron spin coupled with a nuclear bath, 
Eq. (jl5j) . contains two terms. First one is a phase factor due to spin precession in the 
external field plus the Overhauser field. The second term is due to electron-nuclear spin 
entanglement. In ensemble measurements, the inhomogeneous distribution of the nuclear 
Overhauser fields typically leads to a fast ensemble dephasing time Tj*. This complicates 
direct observation of spin decoherence in FID. To remove the undesired phase factor one can 
use, for example, a Hahn spin-echo setup. However, it should be noted that the magnetic tt 
pulse affects the entanglement term also.— All the terms evaluated in the previous section 
can be calculated for the echo setup straightforwardly. For example, the second order term 
with the dipole-dipole interaction. Fig. |2l is 

{Vt) = - Y: {BffnuL^ + 4^^^ 4^^}. (16) 

. 'T"^ , '-^ii ^i-i ^i-i 

«=Tj=i ^ ^ ^ 

Moreover, all the diagrams with the hyperfine-mediated interaction are cancelled for the 

echo. 

As an example we apply the developed technique to a model system, a phosphorous donor 
in a Si crystal, that was first studied long ago.— Recently, interest in it has been renewed by 
a proposal for quantum computation.— The '^^P is a shallow donor with the effective radius 
of the electron wave function R^s ~ 25A;51 Therefore, the bounded electron covers many 
host lattice sites. The nuclear spin bath is represented by a system of randomly distributed 
^^Si isotopes (/ = 1/2). The natural ^^Si isotope concentration is c(^^Si) ~ 4.7%. At 
temperatures ~ IK and magnetic fields ~ 0.1 — IT the major mechanism of spin echo decay 
in this system is a nuclear spin spectral diffusion 

We simulated numerically the processes shown in Figs. I2I9I for a single FID (7r/2 — t— 
measurement) and spin echo setup (7r/2— 1/2 — tt — t/2— measurement). For FID we factorize 
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the shift due to the nuclear Overhauser field out and concentrate only on the decoherence in- 
duced by electron- nuclear dynamical entanglement. Such calculation becomes relevant when 
the inhomogeneous broadening can be filtered out, e.g., by the method discussed in Ref. l58. 
The contact hyperfine constants for the system were approximated using the effective mass 
theory envelop function.^ We also assumed that the phosphorus nuclear spin contributes to 
the frequency shift only, because of large difference in 7— factors of '^^P and ^^Si nuclei. In 
simulations of dipole-dipole contributions we generated an initial nuclear spin configuration 
in a Si lattice within a sphere of radius 5i?efr about the donor. Nuclear spin bath was as- 
sumed un-polarized. Then we selected randomly a spin-up site with its surrounding within a 
sphere of a radius 5a, where a = 5. 43 A is a Si lattice constant. We calculated contributions 
of all possible configurations of a given central spin with its surrounding. This procedure 
was repeated for 10^ times and results were normalized to the total number of spin-up within 
the whole simulated volume. For the hyperfine-mediated interaction we averaged over 10^ 
randomly generated spin configurations within the whole volume. The simulation was done 
for 100 different configurations of nuclear spins to account for initial state dependence of the 
decoherence process. The results were checked for convergence with changing parameters of 
the model. 

For the electron spin free induction decay the real parts of different diagrams. Figs. 12191 
averaged over spatial and spin configurations of ^^Si are shown in Fig. ^| At very short 
times if: < maxjA*^^}^^ the second order terms^ are ~ t"^, and V.^'^ ~ f^. In Fig. [TUT a) 
we show the crossover from the short-time behavior to an intermediate regime with time de- 
pendencies V^^ ~ t^, ¥2'^ ~ t^'^. The dispersions in the exponent at intermediate timescale 
is ~ 5% depending on spatial positions of ^^Si near the P donor and different initial con- 
figurations of nuclear spins. It also include errors due to a finite simulated volume. Unlike 
the electron spin in a quantum dot2£ the short-to-intermediate regime crossover in Si:P is 
more noticeable because of the stronger confinement and inhomogeneity of the electron wave 
function. For external magnetic fields H < 0.1 — IT the hyperfine-mediated term determines 
short-time spin dynamics, see Fig. llOf a). However, it can be efficiently suppressed by in- 
creasing the field. Moreover, it cancels completely in the spin echo. All the dipole-dipole 
fourth order terms develop on the timescale of the order of several msec. Fig. Iior b). How- 
ever, on this timescale the electron spin coherence is completely destroyed by the (V2) term. 
This slow development of high order spin correlations is consistent with an experimental 
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measurement s.~ At a very short times the fourth order terms are ~ (not shown in the 
figures). This dependence changes to V^'^'^ ~ t^'^ at longer timescale. The contribution of 
the four-spin diagram is about an order of magnitude larger than other fourth order terms. 
However, we abstain from attributing it to 1/Z expansion because there is no such difference 
between two and three-spin-fourth order terms. Probably, this is because the studied system 
of nuclear spins is dilute and Z is of order of unity. In Fig. ^2 we show the total second and 
fourth order contributions (see inset) in the exponent for FID and spin echo setups. One can 
see that a vr-pulse reduces electron spin decoherence. This effect is an analog of suppression 
of a spectral diffusion considered in phenomenological models.-^^*^ 

For the spin echo setup we obtain the time dependence of V^'^'^ term comparable to that 
was calculated in Ref. |27l However, the fourth order terms in our model do not show the 
non-monotonic behavior. 

We emphasize that in the Si:P system the non-contact hyperfine coupling between the 
electron and nuclear spins, not considered here, produces noticeable effects on electron spin 
dynamics.— >^^»^ To suppress these effects a high external magnetic field is required.-^ More- 
over, in spin echo measurements on macroscopic samples the dipole-dipole interaction be- 
tween electron spins causes an instantaneous diffusion.^^i^ This effect is beyond the scope 
of this paper on a single electron spin. For comparison with the experimental results given 
in Ref. ^ we account for the instantaneous diffusion with a phenomenological exponential 
decay, e~*/*^^ . In the simulation we take the phenomenological relaxation time for instanta- 
neous diffusion to be t/^ = 1.1 msec, obtained in Ref. Is. The results are shown in Fig. [T21 
Besides the echo modulations due to non-contact dipole-dipole interaction observed in the 
experiments we still have a moderate discrepancy. Because the magnetic field used in the ex- 
periment was of the order of 0.3 T, the contribution of the hyperfine-mediated terms should 
be small on the timescale of the echo decay (0.1-0.5 msec). We attribute this discrepancy to 
the effective mass approximation for the electron wave function. More detailed comparison 
should be done after echo modulation and instantaneous diffusion effects are accounted for 
by the theory. We leave it for further studies. 
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V. CONCLUSION 



We developed a field theoretic approach to evaluate dynamics of an electron spin inter- 
acting with a nuclear spin bath in a high field regime. The approach provides a better 
understanding of the difference between stochastic models of an electron spin spectral dif- 
fusion and dynamic models of spin decoherence in the presence of the nuclear spin bath. It 
also throws light on the problem of reversibility of spin dynamics. The approach is based 
on a conventional diagrammatic technique utilized in the study of Heisenberg ferromagnets. 
The scheme allows for analytical evaluation of different processes contributing to the elec- 
tron spin evolution. We show that electron spin dynamics in a nuclear spin environment 
can be factorized into a free precession in the Overhauser field and more complex dynamics 
due to an electron-nuclear spin entanglement. The latter can be evaluated using a linked 
cluster expansion procedure. The exact analytical expressions for second order and some 
high order processes are given. We show that spin decoherence of a P donor electron in 
a single Si crystal is mostly controlled by nuclear spin pair excitations at sufficiently low 
temperature and high magnetic field. Contributions of high order processes are small and 
can be neglected on a timescale up to several msec. A magnetic vr-pulse flipping electron 
spin slows down the decoherence process. The simulated results are in fairly good agreement 
with experimental measurements of spin echo in macroscopic samples. 



APPENDIX A 

We briefly summarize the spin diagrammatic rules. Unlike Ref. where the formalism 
was used for a statistically mixed state, our brackets correspond to a pure state specified 
by initial conditions. Averaging over a thermal ensemble in our case would mix dynamical 
contributions of electron-nuclear spin entanglement with effects of statistical distribution of 
Overhauser fields. For the sake of simplicity here we assume that / = 1/2. However, the 
approach can be extended to a more general case.— 

The matrix element in Eq. (jH)) can be written as a product of brackets corresponding to 
single sites. At each nuclear spin site j with a given initial state, |j) (= | t) or | |)), an 
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expectation value of time ordered spin operators, 



(j|r{r(ti)/-(t2).../+(4)/^(4+i)...F(t™)}|j) (Ai) 

is zero if numbers of /"'" and I operators are not equal. Otherwise, we evaluate it with the 
Wick's theorem.— Using the spin commutation relations, [J^,/^] = —2/^ and [/^,/+] = 
the bracket ()A1|) is transformed to the form where an operator /+ (J^ would serve as well) 
is in the first position 

{j\f(^^{t)i\t,)\j) = {j\[nt,),i^t)]i^{t,)\j) + (j|/+(t)r(to/"(t2)ij) (A2) 

or in the last position 

(j|r(ti)/^^(t2)|j) = [i^{h),i^{t)]\j) + {j\nt,)i^{h)i^t)\j). (A3) 

Depending on the initial state operator, /"^ is moved to the right if |j) = | t) or to the left 
if |j) = I i)- After applying 

j+(t)| T) = 0, 

^ ^1 (A4) 

(i \i^t) = 0. 

a product of m spin operators is expanded into a sum of products of m — 1 operators. This 
procedure is repeated until only a product of P operators is left. The latter term is evaluated 
directly. As a result, the bracket ()A1|) can be written in terms of all possible contractors of 
operator with /~ and P using 

[J°(ti),/+(t)]=e-(*-*^)[/",/+],„ (A5) 

where the latter commutator is taken at time ti. Unlike contractions of bosons or fermions 
a commutator of and either /~ or P is an operator that can be used in a next pairing. 
For example, in 

f^^'i^rih) (A6) 
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the operator I~^(t), firstly, is paired with P(ti) with resulting I~^(ti) operator paired with 
I^{t2). Another specific example of spin pairing is a locked term.— In 



1 



(A7) 



/+(t)/-(ti)/+(t')/-(t2), 



the operator I^{t) is paired with I~{ti), then I^{t') is paired with the resulting operator 
P{ti) and finally with I^{t2). By the locked term here we mean a term that contains an 
operator I with three contracting lines. The role of such terms in LCE we discuss in 
Appendix |B| 

In diagrams we depict /"'" vertices by points, /~ and P vertices by open circles and 
interaction terms by wavy lines. The Green's function, defined as 



is shown as a line with an arrow propagating from to . An initial spin state determines 
Green's function time evolution. One can see that the Green's function propagates back in 
time (arrow points opposite to the time arrow) if the initial spin state is |, and forward in 
time, if the state is [. Although I~ and P vertices are depicted by same symbols there is 
a topological difference in their appearance in diagrams. P vertex can be either separated 
from any Green's function or connected to one incoming and one outgoing Green function. 
I vertex can have one incoming line or two incoming and one outgoing lines. /"'" vertex 
always has one outgoing Green function line. For interactions, a wavy line connecting two 
circles corresponds to PP term, while a line connecting a point and a circle corresponds to 
I^P term. 

The diagram representation can be easily translated into Green's functions. For example, 
the second order flip-flop term. Fig. |21 can be written as 



The coefficient A in front of the sum accounts for a number of equivalent diagrams. In the 
particular case it is equal one. The coefficient (—2)^ appears from two contractions of 
operators. Analytical expressions for diagrams are dependent on initial spin states. For 




(A8) 
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example, pairs |t or || give zero contribution to (V2) term, see Fig. |21 while contribution 
of ti and It pairs are equal. We usually omit this configuration dependence in graphic 
representation. For example, three spin diagram. Fig. |21^c) corresponds to two possible spin 
clusters t|J, and Iff fFig. I13j) that have different analytical expressions. In general, one can 
distinguish between configurations that can be transformed to each other by changing order 
in spin counting or by rotation or inversion of the coordinate system and configurations that 
are distinct. First type of configurations are t|J. and IfJ,. If we start counting spins from the 
I site these configurations are the same, and they have equal analytical expressions. In the 
case given in Fig. El the two spin configurations are connected by the inversion operation. 
An analytical form of the second diagram at the right hand side can be obtained by changing 
signs at all frequencies uj in an expression for the first diagram. Distinct configurations are, 
for example, titi and ttii contributing to the same forth order ring diagram. Fig. Eff). 

APPENDIX B 

Here we consider in more details a physical origin of locked terms in LCE. For the sake 
of simplicity we assume that the interaction is 

V{t) = Y,B,,lt{t)Ij{t), (Bl) 

and we evaluate 

{n\T{e-'^'^^^''^''''}\n) = {n\l + {-i) f V{t')dt' + {-if /2\ f f T{V{t')V{t")}dt'dt" + ...\n), 

Jo Jo Jo 

(B2) 

where n denotes a spin configuration. The locked diagrams appear in the fourth order 
contribution and correspond to two-spin and three-spin excitations. A two-spin-fourth-order 
correction can be written as 

3(-.)V4!$^4 r r r f T{m{h)i;{t,)it{h)i;{u)\^ (bs) 

JO Jo JO JO 

''J 

{j\Ij{ti)I+{t2)Ij{h)lf{U)\j)}dtidt2dhdU, 
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where we have separated operators corresponding to different spins. The coefficient 3 in front 
of Eq. ()B3|) accounts for possible choices of i and j. Eq. ()B3|) is nonzero only if i =t, j =i 
(for i =l,j =t we just change order in counting of spins and get the same configuration). 
The only possible time ordering in this case is ti > ^2 > ^3 > ^4 or ti > > > t2. The 
integrand for both sets of time ordering is the same and equal e*'^'-'*-*^"*^"'"*^"*''^. The same 
result we can get with the spin diagram technique and diagram equations given in Fig. |H| 
Firstly, we expand the time ordered product of spin operators in Eq. ()B3|) in terms of all 
possible contractors as discussed in Appendix ^ There are two unlinked terms of the form 



/+(ti)/r(t2)/hi3)/r(i4)/7(^i)^;(^2)/7(^3)/;(^4), (b4) 

two linked terms corresponding to the ring diagram, given in Fig. WIS) 



I 



irih)ir{t2)int3WiU)Iiiti)i;ih)Iiits)I^iU), (B5) 
eight single-site locked terms. Fig. Hl^g), of the form 



n Ir 



1 



j+(ti)/r(i2)/+(t3)/-(t4)/-(ti)/;(t2)/-(t3)/;(t4), (b6) 

and four two-sites locked terms. Fig. ID^h) 



Eq. ()B4|1 corresponds to the series expansion of the second order contribution. However, it 
allows some non-physical states, because the contractions in it specify time ordering ti > t2 
and ^3 > ^4 only. By direct evaluation, one can show that the locked terms, Eqs. flB6IB7|) . 
modify it according to the diagram equation given in Fig. IHa) to get a correct time orde ring 



ti > t2 > H > tii. It is a correction of overlapping of spin excitations discussed in Refs 
The locked diagram is the price we pay to obtain LCE for spins. This shows that the 
procedure to get an exponential form of a qubit decoherence is not as transparent as it was 



suggested in Ref. 



27 
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APPENDIX C 



Here we provide explicit analytical expressions for two and three spin diagrams given 
in Fig. El The first term is denoted as V43, the three-spin part of the second term is V^u 
and the two-spin part of the second term plus the third term is V42Z. We do not include the 
coefficients 1/2 indicated in Fig. El The analytical expression for the last, four-spin, diagram 
is lengthy though its evaluation is not complicated. We give terms for the hyperfine-mediated 
and dipole-dipole interactions separately. 

Contributions of the hyperfine-mediated interaction are: 

a) the three spin ring with overlapping corrected, 



b) the diagram compensating a single site overlapping. 



{ym) = 2 2^i^,,i^ifc|2^t— -2- + , ,3 , ,2 , , (<^2) 

I _ ^2i{uJik+'^i3)t 
} + {j ^ k}, 



c) the compensation of a double site overlapping (two diagrams) 



(«> = -E4{-2'*i^^ + ^^^ ^.y (C3) 

The same diagrams for dipole-dipole interaction only have the following expressions: 
a) the three spin ring with overlapping corrected. 



(>/43 ) = - 2^ ^*i^jfci2«H 7-27-7^ -27^ ) + 4 ^,3,, ,,2 (^4) 

^^ij^^ik^^jk ^ij^jk ^ij^ik^jk 

,2 , ,2 ^ , , , , , ,2 , ,3 , ,2 , , , , , ,2 ' , ,2 , ,3 J ' 



Hfc^jfc ^ij^ik^jk ^ij^jk ^ij^ik^jk ^ik' 

b) the diagram compensating a single site overlapping, 
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} + {j ^ k}, 



c) the compensation of a double site overlapping, 

nrddx o4 ro • 2(1 - e*^"^*)2 - 3 e^^-^^^* - Se^*^"^* + 12e2*'^'^* + Se^"^'^* - 13. 
(1/4^2) = -2^^^Ji^^t -3 + TT-i 1- 

In equations {j k} means the same expression with interchanged j and k indexes. 
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FIGURE CAPTIONS: 

Fig. m Schematic representation of two-branch evolution. 

Fig. 121 Second order diagram corresponding to nuclear spin pair flip-flop processes. 
Fig. |3J Set of third order diagrams 
Fig. 1^ Fourth order diagrams. 

Fig. El Second order contribution to decoherence in LCE that includes a series of even 
orders of a conventional perturbation expansion. 

Fig. |BJ Renormalization of a Green's function by /f /J interaction. 

Fig. 13 Compensation of spin excitations overlapping in time by a locked diagram. Time 
overlap in one site (j). 

Fig. IS) Compensation of spin excitations overlapping by a locked diagram. Time overlap 
in two sites (i and j). 

Fig. El The total fourth order contribution of nuclear spin excitations to electron spin 
decoherence. The flrst diagram corresponds the ring propagation of spin excitation in three- 
spin clusters tii and |tt. The second, single-site locked, diagram corrects overlapping of 
excitations. It contributes spin dynamics of three and two-spin clusters as shown in Figs. U\ 
and IHl The third, two-site locked, diagram corrects overlapping of two spin excitations. 
Fig. ISl The last diagram describes four-spin ring excitations in spin clusters ittt) Till) 
nn and TTii. 

Fig. ^1 Contributions of linked clusters to the electron spin decoherence: (a) Crossover 
from the short-time evolution to the intermediate regime, (b) Development of the fourth 
order terms in the intermediate time regime. V2 corresponds to the spin pair flip-flop process 
in Fig. 12 V43 - the three-spin excitation, the flrst term in Fig. V4/1 - the single site 
spin locked diagram, the three-spin part of the second term in Fig. El V412 - correction of 
overlapping of two-spin excitations given in the bracket in Fig. jHl^a), and V4 - the four-spin 
ring diagram, the last term in Fig. IHl 

Fig. ^2 Time dependence of second and fourth (inset) order terms in the exponent. FID 
vs. spin echo decay. 

Fig. ^1 Spin echo decay in comparison with the experimental data.^ Orientation of the 
external magnetic fleld is along (100) direction in the crystallographic axes. 
Fig. ^1 Representation of a third order diagram using cluster conflgurations. 
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